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Abstract 

We derive new bounds on the Dirac mass of the tau and muonic neutrinos. 
By solving the kinetic equation for the rate of energy deposition due to hehcity 
flipping processes and imposing the constraint that the number of effective species 
contributing to the energy density at the time of nucleosynthesis he Ak^ < 0.3, 
we find the bounds m^^ < 150 KeV and m^^ < 190 KeV for Tqcd = 200 MeV. 
The constraint Ak,^ < 0.1 leads to the much stronger bound rrii, < 10 KeV for 
both species of neutrinos. 



If neutrinos are massive and stable, then the mass range 40 eV < m,y < 2 GeV f!], § is 
excluded by measurements of the age of the universe. However, if the neutrinos are unstable 
then these bounds will not apply and we must search elsewhere to find mass restrictions. The 
present laboratory bounds on the muonic and tau neutrino are 160 KeV|p and 31 MeV |Q 
respectively. Better bounds may be had if the lifetime of the massive neutrino is > (9(100) 
sec. In this range of lifetimes the energy density due to the massive species may come 
to increase the expansion rate of the universe and thereby increase the primordial helium 
abundance ^j. In this case the mass ranges 0.5 < -m < 35 MeV and 0.3 < m < 35 MeV were 
found to be excluded for Majorana and Dirac neutrinos respectively [^]. The mass bound for 
Dirac neutrinos was derived assuming that the right handed species is in thermal equilibrium 
below the QCD phase transition, for m > 300 KeV. In this letter we show that even if the 
right handed neutrinos do not enter thermal equilibrium at temperatures below Tqcd, their 
out-of-equilibrium production rate is strong enough to yield a much more stringent lower 
bound on the Dirac neutrino mass. 

The rate of production of right handed neutrinos is roughly given by, r_ oc {m^/ EYG\T^ 
where E is the energy of the neutrino and Tw is the normal weak interaction rate. Thus, 
the "wrong" helicity statesQ are produced more efficiently at higher temperatures. If the 
wrong helicity states decouple from thermal equilibrium above Tqcd, their energy density 
will be greatly diluted in comparison with the coupled species as a consequence of entropy 
conservation. Such neutrinos would appear to be cosmologically safe. A bound based on the 
decoupling argument was first derived in ref. [§], where it was shown that the Dirac mass 
must satisfy m^ < 300 KeV, assuming that the QCD phase transition temperature is 100 
MeV. The result of ref. 0, was refined in ref. 0, where a more careful computation of the 
relevant scattering processes was performed and much weaker bounds were claimed. 

In ref. a different bound, m^^ < 420 KeV was obtained by considering the out-of- 
equilibrium production of wrong helicity states through the pion resonance 77 ^ 7r° ^ ^+^+5 
and imposing the constraint Ak^, < 0.3. Unlike the bound discussed first above, this result 
is insensitive to Tqcd, as will be shown below. 



^In the rest of the paper, when we refer to "wrong" hehcity states, we mean right handed neutrinos and 
left handed anti-neutrinos. 



In this letter we give a more accurate analysis of the production of light wrong helicity 
states in the early universe. We will account for all sources of wrong helicity neutrinos, 
including the population that decoupled above Tqcd as well as those generated from the 
following out-of-equilibrium processes: 

hh -^ hi^fi{T)+ (1) 

Moreover, we will accurately solve the kinetic equation for neutrino energy density generated 
out-of-equilibrium at temperatures below the QCD phase transition. Our bounds will be 
both more precise and much stronger than any found in earlier literature. 

We should mention here that strong bounds on the neutrino mass, ttIj, < 10 — 20 KeV, 
were obtained from the consideration of the cooling of supernova SN-87 [|TU|, |TT|. While 
these bounds are indeed more stringent than the bounds obtained in this paper, they are in 
a sense model dependent. Supernova bounds come from calculating the rate at which energy 
is drained from the supernova assuming that the right handed neutrinos are sterile and will 
free stream from the core. However, since the right handed neutrinos are singlets, it is quite 
possible that they interact in yet unknown ways |]12|.0 Furthermore, with better data it 
may be possible that the nucleosynthesis bounds will eventually become more stringent than 
those obtained from supernova considerations. 

We will first consider helicity flipping scatterings z/^(r)-^ — > i'^(r)+^, and annihilations 
11 -^ t'/x{T)+^'/x(T)+ where / is a lepton (not the tau however) and sub-± refers to the helicity 
of the neutrino. The amplitude for t-channel elastic scattering process is given by 

G 

A{iyl; u+l) = -^iy-fa{^ - i^)p+1{cv - Ca1-,)iJ, (2) 



^ If the right handed species is trapped inside supernova, thus obviating the bound derived in ref . ||l^, ^ , 
it would hkely contribute to the energy density at the time of nucleosynthesis. To avoid this, physics far 
beyond the standard model is needed p3[ . 



where, for example in i/z/-scattering cy = ca = \ and in i/e-scattering ca = 2 sin^ Ow — \ and 
cy = |. The projection onto the wrong hehcity state is accomphshed by using 

^ zy+zy_ = i(^'+m,)(l + 75S), (3) 

spins 

where the four-vector s is given by the expression: 

s = —(k',uj'n'). (4) 

Here the hehcity eigenvalue A = ±1 determines the helicity state of final state neutrino, k' = 
\k'\ is the magnitude of the spatial component of the neutrino momentum, u' = (k''^ + m1Y^'^ 
is its energy, and n is the unit vector in direction of k'. 

The amplitude squared for the production of the neutrinos with the wrong helicity in 
elastic scattering is 

|A(z//;z/+Or = l6Gl{{cv±CA)\p-k){p'-l) 

+ {cv T CAf{p' ■k){p-l)- (4 - cl)mf{k . I)} (5) 

where p, p', are the four-momenta of the initial and final fermions, and k, k' are the momenta 
of the neutrinos i/_ and z/+ respectively. The upper (lower) signs refer to scattering off par- 
ticles (antiparticles) and we defined the four-vector / = {uj' — k'){l, —n) = (m^/2A;')(l, —n). 
The amplitude for annihilation process may be obtained from (^) by simple crossing relation. 
One may see that for the scattering off of neutrinos, jAp goes like the neutrino mass to the 
fourth power in the center of mass frame. It is because of this seemingly slow rate that these 
processes were neglected in ref. [0. However, this suppression is not present in an arbitrary 
frame and generically the production of the wrong helicity neutrinos is suppressed only as 
m^. This is a consequence of the fact that the helicity is not a Lorentz invariant quantity. 
Indeed, given a particle with definite negative helicity, a boost in a direction orthogonal to 
its direction of motion with velocity v will generate an admixture of the positive helicity 
equal to m'^v'^/AE'^. 

The kinetic equation for the energy density of the right-handed neutrinos has the form 

dpu ,rT ^ 1 f (fp(fkd^p' d^k' ,^4, , , ,,x 

-^ + 4:Hp^ = Ct = - — rz —^ uj'6^(p + k-p' -k') X 

dt ^ (27r)8 J 2E 2u 2E' 2uj' ^^ ^ ' 

x{|Ar/.H/K^)[i-/z(^')] + •••}' (6) 



where the elhpses refer to the remaining elastic scattering (off antiparticles) and annihilation 
channels as well as to the contributions from the pion resonances which will be treated 
below. Note that we omitted the Pauli blocking factor for wrong helicity states, which are 
not in thermal equilibrium. We also neglected the inverse reaction of right-handed neutrinos 
producing the normal ones. This is justified because the density of the wrong neutrinos has 
to be small. For the same reason we can compute all contributions to p^, separately and add 
the different contributions in the end. 

If we assume that the phase space distribution functions are of the equilibrium Boltzmann 
form, / = exp{—E/T), we can compute the collision integral on the r.h.s. of equation (^) 
analytically in the massless limit. For example, the matrix element for a t— channel reaction 
in (|^) gives 

Ctimj = 0)mb = ^\ff\ icv + CA? + l{cv - ca?) (7) 

We have also computed the collision integrals numerically retaining the effect Fermi statistics, 
and found a 21% suppression relative to the case of Boltzmann statistics. Furthermore, it 
was found that this suppression is nearly the same for all the reactions. 

We took into account the following scattering processes: t'^-Z/g — ^ t'^+^^e; ^^i-^fi -^ 
v^+v^, v^-^r -^ T-'n+Vr, v^^c" -^ t'lj+c", v^-P~ ^ z^^j^/i^, with the corresponding channels 
for scattering off antiparticles and the annihilation channels. Similar channels exist also 
for z/^4.-production. However, the reaction channels involving muons are different, due to 
the contribution from charged current scattering. Moreover, in the muonic case there are 
three more pure charged current channels that do not exist for z/T-+-production: p~V(. —>■ 
u^+e' , fi~e~^ — > t'^+z/g and e'^i'e -^ ^n+l^^ ■ There is also a contribution from scattering off 
of pions, which we found to be small and will be neglected. 

Equation (J^), modified to include Fermi statistics, can be used for all relevant reaction 
channels except for those involving muons. For the muonic channels we computed the colli- 
sion integrals numerically as a function of temperature. Altogether we can write the kinetic 
equation as 



27r5 
where the Fermi correction factor cp ^ 0.79 and the effective number of reactions A^eff 



p, + AHp, = ilfl^cpAT^ffT^ (8) 



(normalized such that t-channel scattering off neutrinos gives a contribution 1), for example 
in z/^-case is given by 
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iVcfr,. = (y)-{l+ E {c'vf + c%) + {4^ + cl^ + 2)d,{T)}. (9) 

The first '1' comes from scattering off of muonic neutrinos and cy^ = 2 sin^ 9w + |, c^^ = | 
(see also below equation (|^)). Using sin^ 6'^^ = 0.23, we find 

A^eff.M ^ 3.92(1 + 0.81t/^(T)) 

A^eff.r ^ 3.92(1 + 0.06rf^(T)), (10) 

where the fit functions di{T) are normalized to become unity for very large T. The sup- 
pression due to the muon mass becomes significant rather late: for T = m^/3 we still get 
d^ ~ 0.42 and dr — 0.35. Hence the muonic interactions are important in deriving the mass 
bounds. For T = m-^/lO however, we find t/^ ~ 4.7 x 10^^, which clearly validates our neglect 
of reactions involving rs. 

If all the particles in the plasma were massless we could easily integrate the equations 
(§-[Ty) since in this case the temperature and the time are related as T^t = const. However, 
since the contribution of massive pions and muons to the energy density is essential, we 
should proceed more carefully. If we neglect the tiny contribution to the pressure density p 
coming from the wrong helicity states, the covariant energy conservation law p = —3H{p+p) 
reduces to the usual conservation of entropy of interacting species, and we get the standard 
relation between time and temperature |T| 



dt _ .45 .i/a ^p; 1 / T dhA 

dT Mvrs^ r3 gl/^^T) [ SKdTj- ^ ' 

Here g^, and h^, are the effective numbers of energy and entropy degrees of freedom, defined 
through the total energy density p = ^g^(T)T'^ and the entropy density of interacting species 
sj ^ ^-iK{T)T\ 

Equations (|8|-pTD are easily solved to yield a simple integral expression for the relative 
energy density ri = Pu+/pu- = Pu+/{77t'^T'^/2A0): 

/ m^ \2 /-T^QCD / 10.75 X^^V 10.75 X^^V T dhA , , , ,, 



Sub-1 means that we take into account only weak scatterings; pion decays are treated sep- 
arately below. Tqcd is the temperature of the QCD phase transition and the coefficients Oj 
can be read off the equation ([lOD , a^ = 0.81 and a,- = 0.06. Equation (|l^) incorporates the 
effect of the changing energy density due to annihilations of heavy particles, as well as the 
statement of entropy conservation. We have integrated (|12D numerically and found that (for 
^QCD ~ 100) to a very good accuracy r^s follow the simple linear fits 

ri^ = (0.53 + 2.25 T^^^^) (m,/MeV) 

rir = (0.08 + 1.34 T(J°cd) (^./MeV), (13) 

where TA^q is the QCD phase transition temperature in units 100 MeV. A word of caution 
is in order here: In equation (|12D we assumed that g^ and h^, reach a constant value 17.25 
well above pion mass. This indeed seems to be the correct choice |T^, but should they 



instead increase significantly at temperatures between m-Tr and Tqcd, then equations (|T3D 
would slightly overestimate ri for Tqcd ~ 0(200) MeV. 

The decays tt^ — > ui? and tt"^ -^ ^v^ also give significant contributions to the energy 
density of right-handed neutrinos. The decay of vr*^ produces both v^ and Vr while vr^ 
produces only z/^. The former was considered in ref. while the latter, to best of our 
knowledge, has not been accounted for to date. Both of these processes have the property 
that their contribution to the energy density is insensitive to the temperature of the QCD 
transition because the dominant contribution to production of the right-handed z/ occurs at 
T K, niT^/b. Let us ffist consider the decay vr*^ -^ ^+^+- This width was calculated in ref. 



I|, |T§ as 

r(7rO ^ z/+zy+) = GlflomlmjM, (14) 

where Jt^o c^ 93 MeV. When computing the 7r° collision term, we will assume the Maxwell- 
Boltzmann statistics. This is an excellent approximation here and moreover we are able 
to compute the collision term analytically. We find that the energy density of neutrinos 
produced in 7r°-decays satisfies the equation 

p. + AHp^ = ^^f^,T'xtK,{xo), (15) 



where xq = r/i^o/T and K2{x) is the usual modified Bessel function. This equation can be 
integrated analogously to the scattering contribution (cf. equation (p!2D) with the result 



r2 = 1.21(m^/MeV)^ (16) 

which to the given accuracy is independent of Tqcd- Quite similarly we find for the decay 

P. + ^Hp, = ^|^/'±^'4(i - 2/'){(i + y')K2{x^) - y'Koix+)}, (17) 



where /,r± — 128.4, a;+ = mT,±/T and y = m^/mjr±. When integrated, (p!?!) gives 

rs = 1.15(m^/MeV)l 



At temperatures above the QCD phase transition, the plasma consisted of free quarks 
and gluons in addition to leptons and photons and the effective number of degrees of freedom 
g^ was slightly below 60. Since the efficiency of production of right-handed neutrinos is larger 
at higher temperatures (at least up to the temperatures of the order of the W and Z boson 
masses) we expect them to be in thermal equilibrium if their mass is above (9(10) KeV. The 
energy density of these neutrinos will be diluted in the course of Q CD-phase transition and 
due to annihilation of massive states in the plasma. Using the entropy conservation we may 
estimate their relative energy density as 

rQCD = (10.75/60)^/3-0.10. (19) 

There is one additional effect that needs to be taken into account. Namely, the fact 
that energy density of a massive species will scale differently than that of a massless species. 
Thus, the there will be additional contributions to the relative energy density coming from 
population of left handed neutrinos as well as the small population of right handed neutrinos. 
Therefore, the total relative energy density may be written as 

rtot,u = f{m,) + (Y. ^i,- + 0-1) (1 + /(^-))- (20) 

i 

The function fim^) was found using the results from paper |^, and is given by 

/(m^) = 1.31m^ -f AMml - A.blml = 1.31ml + Sf, (21) 
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and rriy is given in units of MeV. In this way we find the bounds 



/^^max_Q_^0-5 



1/2 



m^ < '■ — Tan^ MeV, (22) 

1^ 4.33 + 2.25r^0c°D / ^ ^ 



/^^max_Q^^0-(5^~ 



1/2 



where Sf, = 1.105/ + (2.89 + 2.25T^^^)f and 5^ = l.lOSf + (1.29 + 1.34T(^°cd)/- Due to 

smallness of 6i ~ 0{ml), the constraint equations (|2^) and (0) are easily solved for any 

values of AA;,^ and Tqcd- For example, imposing the constraint /S.ky < 0.3 [jl^, leads to the 

bounds 

< r 170 KeV Tqcd = 100 MeV 



170 KeV 
150 KeV 


Tqcd 
Tqcd 


= 100 MeV 
= 200 MeV 


210 KeV 
190 KeV 


Tqcd 
Tncn 


= 100 MeV 

= 200 MeV. 



"^ ~ I 190 KeV Tqcd = 200 MeV. ^ ^ 

Let us finally mention that there were claims recently that the effective number of extra 



neutrino species is not larger than 0.1 ||18[. If this is indeed the case the right-handed 
neutrinos should decouple before or near electroweak phase transition. Using this bound the 
limit on both the i>r and u^ masses is, conservatively, 10 KeV. 

A.D.Dolgov is grateful for hospitality to the Department of Physics of University of 
Michigan where this work started and to the Astrophysical Group at Fermilab where it was 
finished. K.K. wishes to thank the Finnish Academy for financial support. This research is 
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